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Abstract 

Acker et.al. have recently proposed an economical solution to the solar and atmospheric 
neutrino deficits, in which both are explained by large-angle v e — oscillations, sup- 
plemented by v e decays. We show how to embed their phenomeno logical model into an 
electroweak framework in which global electron and muon numbers, (U(l) e xU(l) M ), spon- 
taneously break at a scale of 1 keV. Despite such a low scale, our model is technically 
natural. The naturalness requirement, together with nucleosynthesis constraints, point to 
the existence of relatively light, largely sterile neutrinos with masses in the MeV range. 
We find a number of potentially interesting experimental implications of these models, one 
of which is an explanation of the excess events that have been found near the endpoints in 
the double beta decay of several elements. One formulation of our model involves a novel 
realization of supersymmetry, for which the new light particles and their superpartners are 
split by very small amounts in comparison with the weak scale. 
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1. Introduction 



The various experimental neutrino anomalies that have come to light over the past 
years fall into two categories according to whether they can be explained simply by a 
particular pattern of neutrino masses and mixings, or whether they require in addition the 
inclusion of new, light scalars at very low energies. The solar neutrino problem [1], [2] and 
the atmospheric neutrino deficit [3] are in the first of these categories, while the observed 
excess of electrons in the double-beta {(3(3) decay of 76 Ge, 100 Mo, 82 Se, and 150 Nd [4] fall 
into the latter group. This basic difference makes it difficult to understand all three as 
being different manifestations of the same type of new physics. 

Recently Acker et.al. [5] have proposed an economical model within which both the 
solar and atmospheric neutrino deficits are simultaneously explained using only electron 
and muon neutrinos. In their approach, the atmospheric v li jv e ratio is depleted by v e — 
oscillations with Am 2 ~ 10~ 2 eV 2 , and sin 2 (26> eM ) ~ 0.8. This region of parameter space 
is unconstrained by the 1MB limits on upward-going muons [6], once matter-effects are 
included for the propagation of z/ e 's through the earth [7]. Besides being depleted by these 
same oscillations, solar neutrinos are further suppressed in this model by neutrino decays 
en-route from the sun: Vh — > V£ + cp. In these decays <p is a scalar particle - perhaps a 
majoron - and the heavier mass eigenstate, Vh, is dominantly v e . It is this additional decay 
that permits a solar-flux signal at Homestake that is less than half of what the standard 
solar model predicts, a suppression that would be impossible with only large-angle two- 
flavour oscillations. 

The one experimental drawback of the model is the energy-dependence that it predicts 
for the solar-neutrino reduction, since neutrino decays preferentially suppress the flux of 
lower-energy p-p-cycle neutrinos that are seen in the Ga experiments in comparison with 
the higher-energy 8 B neutrinos which dominate the signal in CI [8]. This suppression 
arises because the perceived lifetime of the more energetic particles is longer (due to time 
dilation) in the rest frame of the solar system. This prediction is ruled out up to the 95% 
C.L. by the observations, which indicate that the Homestake solar-neutrino suppression is 
larger than that seen in GALLEX. 

The advantage of this picture, on the other hand, is that the use of neutrino decays 
into scalars casts the solution to the atmospheric and solar problems in a way which is 
very similar to the interpretation of the excess (3(3 events as the signal of scalar particle 
emission [9], [10]. This raises the possibility of accounting for all three of these experimental 
anomalies using the same kind of new physics. With this in mind, we choose here to 
reserve judgement on the potential discrepancy with the relative size of the GALLEX and 
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Homestake experiments, in order to further explore whether the (3(3 anomaly can also be 
described. In any case, given the difficulty of these experiments, with the potential for 
unknown systematic errors, we regard this discrepancy as being merely preliminary. 

Our second goal is to come to grips with the fine-tuning problem that plagues all such 
models which involve very light scalars that are relatively strongly-coupled to neutrinos 
[9]. In some models this fine tuning arises because a scalar mass must be directly tuned 
to be extremely small (tens of keV or less). Alternatively, if the light scalar is a (pseudo-) 
Nambu-Goldstone boson, its coupling to neutrinos is typically of order g v ~ m„//, where 
m v is a relevant neutrino mass and / is the scale of symmetry breaking. The requirement 
that g v be large enough to reproduce the anomalous experiments together with the upper 
limit on m v typically implies a small symmetry-breaking scale, / (also tens of keV or less). 
In either case a very small scalar symmetry-breaking scale is required, whose stability 
under renormalization must be established. 

We address these issues by embedding the Acker et.al. theory into a renormalizable 
model of physics at the weak scale and below. Within this framework we are able to find the 
conditions under which the model can solve the atomospheric and solar neutrino problems, 
as well as the (3(3 excess, in a technically natural way. By technical naturalness we mean 
that renormalization corrections to small dimensionless parameters are automatically of 
the same size as, or smaller than the tiny parameter itself. This condition is satisfied by 
all known heirarchies of scale in physics, and so is the bare minimum that might be asked 
of a reasonable theory. It is a weaker condition than 't Hooft naturalness, which would 
also demand that additional symmetries arise as the value of the small parameter goes to 
zero. 

We find that the conditions for naturalness are very restrictive, particularly when 
taken together with what is needed to describe the (3(3 decay excess and nucleosynthesis 
constraints. As we show in detail below, a renormalizable description of the new physics 
requires the existence of new degrees of freedom, which we take to be sterile neutrinos. 
Primordial nucleosynthesis and (3(3 decay argue that these new states should have masses 
larger than of order MeV in scale. This constraint, which tends to drive the new physics 
to higher energy scales, contrasts with the naturalness constraint which prefers the new 
particles to be lighter than 1 MeV. These bounds work against each other, and are con- 
tradictory in the generic case. We find that all constraints can be accomodated, however, 
if the physics of the weak scale should be supersymmetric. Supersymmetry (SUSY) helps 
because it weakens the strength of the naturalness constraints, permitting them to coexist 
with the cosmological and (3(3 requirements. 

We present a renormalizable realization of both the supersymmetric and nonsuper- 
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symmetric scearios here, and find that, for both of these, the combination of all of the 
bounds points to the existence of sterile neutrinos in both the few MeV and the several- 
hundred MeV ranges. The sever al-MeV neutrino is predominantly mixed with z^, and the 
several- hundred-MeV neutrino mixes most significantly with v e . Both types of mixings 
may be amenable to detection through closer scrutiny of the 7r — > \iv decay spectrum, 
as well as in weak-universality measurements of the electron charged- current interaction 
strength. 

2. The Atmospheric and Solar Neutrino Deficits 

We start by reproducing the model of Acker et.al. [5] as the following low energy 
effective theory. In addition to the standard-model (SM) particles and gauge symmetries 
we require the conservation of a global U(l) e xU(l) M electron- and muon-number symmetry, 
which will eventually be spontaneously broken by two additional electroweak-singlet scalars 
which transform under the lepton symmetries in the following way: 

e ~(-2,O) eM ~(-l,-l) (1) 

Two scalars are required here in order to obtain neutrino decays, since a single scalar 
produces Yukawa couplings that are automatically diagonal in a basis of neutrino mass 
eigenstates. 

Under this symmetry only the following SM particles carry nontrivial charge 

e«~(l,0) 

//*~(0,1) (2) 

With these charge assignments, and denoting the SM Higgs by H, the lowest dimension 
terms that lead to neutrino masses and mixing arise at dimension six: 

^2 {\ e ^ e HH(j) e + 2A M Ve^ M M<M • (3) 

This effective description contains all the information that is relevant for atmospheric- and 
solar-neutrino physics, (and as we will see in section 3, the double-beta decay anomaly as 
well). 



(1,0) 

(o,i) 
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Replacing H, <f> e and e/i by their respective f.e.f.'s, v =174 GeV, w e , and gives 
rise to the following mass matrix and Yukawa couplings: 

/ m ee m eu \ ( m ee /w e m eu /w u \ 

m= m ' g= m lw • (4) 

where m ee = X e v 2 w e /A 2 and m eM = \^v 2 w^/ A? . The mass eigenstates and their masses 
work out to be: 

(.) 

where we define c = cos 6> eAt and s = sin 6> eM and 



sin (20 e/ J = — — , r^. 6 

(1 + m e m h ) 2 



The choice sin 2 (26> eM ) ss 0.8 corresponds to m ee pz ±m eM , and Am 2 ~ 10 2 eV 2 then gives 
nih ~ 0.1 eV and me ~ 0.04 eV. (Here Am 2 = m| - mf.) 

The spectrum of light scalars contains two massive states with masses < Wi, as well as 
two massless Nambu-Goldstone bosons 1 which we take to be the (appropriately normalized) 
phases of the two scalar fields: Xe and % eM . Anticipating the result Wi ~ 1 keV that is 
required to produce the desired neutrino lifetime implies that only the massless states are 
relevant to neutrino decays. Their low-energy neutrino couplings are: 



^el 75^e) + ^7 75^ 2 



(7) 



which produces the following total rate for decay, in a frame for which has energy 
E: 

sin 2 (2fl e ,) (Am 2 ) 3 ( 1_ J_\ 
128tt£ m 2 V^e ' 

It is easy to check that this effective lagrangian solves the atmospheric and solar 
neutrino problems in the way envisioned by Acker et.al, as long as the v.e.v. of the scalars 



1 Any masses these states may obtain due to the electroweak anomalies in the global symmetry group 
are completely negligible for our purposes. 
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is about 1 keV, which with the conditions below eq. (6) implies the coupling m e i/wi = 

xy/A 2 ~ 10- 4 . 

3. Double Beta Decay 

The remarkable fact is that these same values that resolve the solar- and atmospheric- 
neutrino problems are also just what is needed to account for the excess high-energy 
electrons near the end point of the double beta decay spectrum. We establish this point 
in the present section. 

The differential decay rate for double-beta decay of a nucleus of charge Z and mass 
number A accompanied by the emission of a light scalar particle is given by [9] : 

dV 

- TF r\yV\ 2 (Q-e 1 -e 2 )[e lPl F(e 1 )}[e 2 p 2 F(e 2 )). (9) 



de 1 de 2 32tt 5 ' 

where G F is Fermi's constant; Q = M(Z, A) — M(Z + 2, A) — 2m e is the kinetic energy 
- typically several MeV — that is available to the final-state leptons; Si and pi are the 
energy and momentum of each of the final two electrons; and F(e) is the Fermi function 
which describes the distortion of the electron spectrum due to the nuclear charge. 

W represents the following integral: 



ij 



d q 

4 



m Ui a ij m Vj - q 2 b 



l 3 



(q 2 + m 2 . — ie) (q 2 + m 2 . — ie) 



w(q°,q 2 ). (10) 



The sum here is over all neutrino species, with V e i representing the strength (normalized 
so that V e i = 5 e i in the standard model) of the e — Vi charged- current interaction, and 
bij are coupling matrices defined by the following Yukawa interaction: C = (bijtfi + 

a ij4>*)lLi / j +c.c. w represents a particular combination of Lorentz-invariant form factors, 
as defined in Ref. [9], which describe the nuclear matrix element of two hadronic charged 
currents. 

For the purposes of analyzing our model, we make the following three simplifying 
assumptions: 

• 1: We neglect the electron mass, m e ~ 0, in performing all phase space integrals, 

• 2: We neglect the Coulomb-distortion factor, F(e) ~ 1, and 
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• 3: We parameterize the nuclear form factors by approximating them by step functions 
in energy and momentum: w ~ Wo <d(q° — E F ) Q{\q\ —p F )- Here p F represents the nucleon 
Fermi momentum and E F = p F /2m N is the corresponding Fermi energy. We fix the 
normalization, wo ~ 4 MeV -1 , by comparing to the observed fifav half lives, and we fix 
p F ~ 60 MeV (and so E F ~ 2 MeV) by requiring that the present upper limit on the half 
life for neutrinoless decay [11], Ti(j3f3o v ) > 2-10 24 yr, imply the corresponding upper bound 
on the v e majorona mass: m Ve < 1 eV [12]. These approximations produce similar results 
to the more detailed calculations of Ref. [9], and are sufficiently accurate to determine 
whether our model can account for the observed anomaly. 

With these choices, and neglecting the neutrino masses in comparison with the nuclear 
scales that are involved, the expression for the matrix element, W, simplifies, in our model, 
to 

(2Try w e 

Using this expression, we tabulate in Table I the values that are required for m ee /w e to 
reproduce the observed anomalous signal. Following Ref. [9] , we integrate the anomalous 
rate only above an energy, E t h, which represents the point above which the contribution 
of the conventional two-neutrino decay, fifizv, is negligible. 



Element Q (MeV) E th (MeV) T? /2 (yr) m ee /w e 

76 Ge 2.0396 1.5 5 x 10 22 8 x 10" 5 

82 Se 2.995 2.2 5 x 10 21 7 x 10" 5 

100 Mo 3.0328 1.9 3 x 10 20 2 x 10" 4 

150 Nd 3.3671 2.2 3 x 10 19 5 x 10" 4 



Table I 

The Effective Yukawa coupling strength that is required to reproduce the anomalous events in double beta 
decay. T^ 2 is the half-life of the anomalous events only, defined to be the total number of events above a 
threshhold value, E t h, for the sum of electron energies, above which essentially only excess events appear. 
m ee /w e is the coupling (defined in eq. (4)) that is required to explain the excess rate. 

It is clear from the table that the required values for m ee /w e are quite consistent with 
the requirements — m ee ~ 0.1 eV and w e ~ 1 keV — that were found earlier as being 
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required to solve the solar- and atmospheric-neutrino deficits. 

Notice also that all four scalar states, and not just the two Nambu-Goldstone bosons, 
contribute to this decay since the energies available, Q ~ MeV, are much larger than the 
typical scalar masses, Wi ~ 1 keV. It is also important for the above analysis that the scale, 
A, of the new physics which is responsible for the effective interactions we are using, c.f. 
eq. (3), is heavier than the nuclear physics scale, p F ~ 60 MeV, of the virtual neutrinos in 
the scalar-emitting (3(3 decay. This is not just a technical requirement because, as we shall 
see, once a renormalizable model is constructed, the (3(3 decay rate becomes suppressed 
if all neutrino states should be light compared to the MeV scale [9]. As a result, an 
explanation of the electron excess within the present framework ultimately requires some 
new degrees of freedom that are not light on the scales that are relevant for double beta 
decay. 

4. A Renormalizable Model 

While the heavy mass scale, A, in our effective theory is large compared to the scale of 
neutrino physics, it is not necessarily large compared to the weak scale. In order to verify 
that our model does not contradict standard electroweak phenomenology and standard 
cosmology (in particular nucleosynthesis) we need to construct a renormalizable model 
for the interactions between the new particles and SM particles at higher energies. This 
is the purpose of the present section. It is also only once this underlying physics has 
been modelled that we may investigate the naturalness of this scheme. We emphasize, 
however, that it is the effective theory presented in section 2 that explains the anomalous 
experimental results and any renormalizable model that leads to it at low energies will do. 

Consider, in addition to the scalar fields discussed to this point, four left-handed 
electroweak-singlet Weyl fields, sf and s^, which respectively carry the charges (±1,0) 
and (0, ±1) under 11(1)6x11(1)^. The most general renormalizable couplings between these 
and the usual standard-model particles are: 

g e s~ip e H + g^s'ijj^H + yS+s+(/> e + /i M s+sj0 eAt + yS~s~</>* + h^s~s~(pl^ ^ 
+ m Se s+s- + m s ^s+s~ + h.c, 

which reduces in unitary gauge to: 

h e h e — 

g e s~u e (v + p)+ g^s'u^v + p)+ ys+s+0 e + h^sts+(j) e ^ + yS~s~</>* + h^s~s~(j)l^ 

+ m Se sts~ + m s ^s+s- + h.c, 

(13) 
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if p denotes the physical SM Higgs. 

Since the lepton-breaking i\e.i\'s, Wi, are so small, we may as a first approximation ne- 
glect them. Then the neutrino mass matrix breaks up into two three-by-three submatrices 
for both the electron- and muon-neutrino sectors. These submatrices have the form: 
















9iV 





9iV 






(14) 



where i = e, fi. Each has one vanishing eigenvalue, and a pair of nonzero ones: ±M S . = 
±^/m| + (giv) 2 . The corresponding eigenvectors are 




and 




(15) 



This gives the following relation between the weak-interaction and mass eigenstates: 

m Si m Si -V2giV 
M Si -M s% 
giv V2m Si 




V2M S 




(16) 



If the heavy Dirac fields, S i , are integrated out, then the effective interaction of eq. (3) 
is obtained with: 

A e gih e 



A 2 M 2 

* e 



and V=WW, 



(17) 



Notice that the more weakly coupled these sterile neutrinos are, the lighter they must be 
to preserve the desired size of m ee and m eM , (recall m ee /w e = X e v 2 /A 2 and m efI /w^ = 

It is worth recalling, in this regard, that if the electron-type sterile neutrino, Sf, 
is much lighter than several MeV, then the predicted anomalous (3(3 decay rate becomes 
suppressed, since in this limit the result is proportional to the v e — v e element of the full 
renormalizable mass matrix, which is zero. This is a special case of the more general 
suppression, discussed in Ref. [9] , of the anomalous double beta decay rate in the low-mass 
limit. We are therefore led to prefer M Se > 60 MeV. 

A useful quantity for comparison with phenomenological constraints is the mixing 
angle which controls the strength of the participation of the sterile states in the charged- 
current weak interactions. From eq. (16) we see that the sterile mixing into the electron 
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and muon charged currents are given by: 

V ^ =v ^ = ^k: and v ^ =v ^ = Tmr; (18) 

These combinations are useful because they are typically constrained to be small, with the 
strength of the bound depending on the mass of the heavy neutrino state, and they are 
bounded from below in the models we consider. Taken together these exclude some of the 
potential mass ranges for the sterile neutrinos. 

In the MeV mass range the experimental constraints come from 7r and K meson 
decays. For example, measurements of r(7r — > eu)/F(n — > \iv) strongly bound the mixing 
into the electron charged current for neutrinos in the particularly interesting mass range 
between 1 and 100 MeV. Such a mixing with a 1 MeV neutrino must satisfy | V e s| 2 < 10 -3 
at the 90% C.L., whereas the same bound for a 10 MeV and a 50 MeV neutrino is 10 -5 
and 5 x 10 -7 , respectively. [13]. Searches for a nonstandard component to K — > ev [14] 
extend this limit up to neutrino masses of 350 MeV. 

Somewhat weaker constraints on sterile neutrino mixing with come from measure- 
ments of the muon spectrum in n — > \xv [15]. For example a 1 MeV neutrino must satisfy 
|^Ats| 2 < 10 -2 at the 90% C.L., and mixing with a 10 MeV neutrino must be less than 
10 -5 [16]. It is noteworthy that beam-dump bounds do not apply to our model, since these 
experiments typically search for neutrino decays into charged particles, and our neutrino 
eigenstates predominantly decay invisibly into neutrinos and light scalars. 

In our model, a lower bound on these parameters starts from the following relation (see 
eqs. (17), (18)): 

\V e sf = W eS +f + \V eS7 f=^. (19) 



e"e 



This, when taken together with the value m ee /w e ~ 10 4 , and the perturbative limit, 
h e < 4tt , leads to the inequality |V e sJ 2 >8x 10" 6 . Thus M Se < 10 MeV. Since these 
searches extend to masses of order 350 MeV, we've excluded M Se in the range 10-350 MeV. 
Likewise |KsU|V^?J = m eM /(«; M /i M ) > 8 x 10 -6 . 

In summary these bounds tell us is that the phenomenologically acceptible parameter 
range is the one for which giV <C M Si ~ rn Si . With this model in hand, we may now 
confront the remaining phenomenological constraints, as well as the naturalness issue. 

5. Cosmological Bounds 

The standard Big-Bang model of cosmology and primordial nucleosynthesis furnishes 
strong constraints on any model which contains light neutrinos and scalars, as do our 
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models [17]. The success of Big-Bang nucleosynthesis sets an upper bound to the number 
of gravitating degrees of freedom at the time when the photon temperature is of order 
T 7 ~ (0.1 — 2) MeV. The more degrees of freedom there are at this point, the faster the 
universe expands, and so the more neutrons are available to be cooked into 4 He. This 
raises the predicted primordial mass fraction, Yp, of 4 He. Since present limits require 
0.22 < Y P < 0.24 at the 95% C.L. [18], and since this agrees with what is expected for 
the standard model particle content, there is little room for new particles to be abundant 
at this temperature. It is conventional to express the resulting bound as a limit on the 
number of neutrinos, 1.3 < iV e fj = N u + 5N < 3.2, where the standard-model value is 
N u = 3, and the observationally preferred value is N v = 2.2. 

Every new spin-half or spinless relativistic particle that is present when the photon 
temperature is T 7 contributes to 5N an amount: 

5N= m/T 7 ) 4 +l^ (T,/T 7 ) 4 , (20) 

fermions bosons 

where Tj is the temperature of particle type l i\ Since Tj = T 7 for particles that are in 
equilibrium with photons, any such particles in addition to the usual three neutrinos are 
basically ruled out. For two light complex scalars, for example, 5N = 16/7, which is clearly 
too large. 

There are two standard ways to avoid this bound, (i) The above expression does not 
apply to particles which are nonrelativistic and still in equilibrium, since for temperatures 
below their rest masses the abundance of such particles is suppressed by the Boltzmann 
factor, e~ m l T . (ii) Alternatively, if a particle decouples from the photon bath at sufficiently 
early times - in practice this means before the QCD transition at T QCD ~ (0.1 — 1) GeV - 
they can be diluted by subsequent reheating of the photons. In this case the ratio Tj/T 7 
can become sufficiently small to suppress this particle's contribution to 5N. 

Only the first of these mechanisms is available to us here, and this only for the poten- 
tially heavy particles such as Sf. Scenario (ii) can not occur in our model since all of the 
heavy and light neutrinos are in equilibrium with one another due to the exchange of the 
light scalars, fa. (In fact for temperatures above T v ~ 2.3 MeV, all of these particles are 
also in equilibrium with ordinary matter because of the weak interactions of the ordinary 
neutrinos.) 

As a result, any heavy particle in the neutrino sector whose mass satisfies m > T v 
satisfies the conditions of loophole (i) above. As T falls below its mass, it annihilates out, 
reheating both the neutrino and ordinary sectors, and therefore not contributing to SN. 
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Any of our new particles that are heavier than a few MeV are therefore cosmologically 
benign, due to their equilibrating interactions involving the light scalars. 

Neither option can apply to the light scalars themselves, however, and these parti- 
cles already contribute unacceptably to SN. Another mechanism is required in order to 
suppress these. One possibility arises if one of the heavier states should have a mass, M, 
that lies below T u ~ 2.3 MeV, but above the temperature, T n / p ~ 0.7 MeV, at which the 
neutron-proton ratio freezes out. In this case, when T falls below M, annihilation into light 
scalars raises the temperature of the neutrino sector, thereby increasing, in particular, the 
number of z/ e 's. As is pointed out in Ref. [19], however, an enhanced v e density at this 
time suppresses the neutron abundance, and so decreases the amount of 4 He that is ulti- 
mately produced. Quantitatively, if n/no and E/Eq respectively denote the ratio of the v e 
number densities and energies before and after annihilation, then the excess z/ e 's (and z7 e 's) 
suppress Yp by an amount equivalent to SN = — 4.65n, where Sn = (n/n )(E / E ) 2 — 1. 

We next attempt to estimate the implications of this effect for our models. Suppose, 
therefore, that N F species of Weyl fermions and N$ real scalars should annihilate when 
T n / p < T < T v , leaving n F and ns fermions and scalars at low temperatures. Then 
the increase in u e energy and density due to the annihilation of these particles may be 
estimated to be (n/n ) = (T a{ter /T heiore ) 3 and (E/E ) = (T after /T before ), where entropy 
conservation requires (T a f ter /T b efore) 3 = [(n F + N F )+j (n s + N s )]/[n F +| n s ]. As a result 
we find 



Sn = 



(n F + N F )+l(n s + N s y b 1 



n F +^n s 



~ 1- (21) 



For example, for the models under consideration there are three light neutrinos and 
two complex scalars at low energies, so n F = 3 and ns = 4. If one of our massive Dirac 
neutrinos, Sf, should annihilate out at this strategic time then N F = 2 and Ns = 0. Then 
this leads to Sn = 0.71 and AT eff = 2.0. This is well within the experimentally permitted 
range, which actually prefers fewer than three neutrinos! 

For future reference we also comment on the consistency of the supersymmetric model 
with cosmological constraints. In such a model the two light scalars and 1.5 MeV Dirac 
neutrino as well as their (approximately degenerate) superpartners must be added to the 
standard three neutrinos. In this case we have n F = 5, ns = 4, N F = 2 and Ns = 4. We 
thus find Sn = 1.2 and N e g = 1.9. There is also the gravitino in the supergravity version 
of the model. It has mass of order 1 keV or smaller and thus satisfies the cosmology bound 
necessary to prevent it from closing the universe. This bound states that TOg rav iti no < 1 
keV or >10 TeV, [20], [21]. 
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In summary, we see that if this is the mechanism for evading the nucleosynthesis 
bound, then cosmology argues for having one of the sterile Dirac neutrinos with a mass 
of around 1.5 MeV, with the second sterile Dirac state being heavier. This heavier state 
should be the electron- type state, 5^, if the double beta decay anomaly is to be repro- 
duced. 5^ must be very heavy in order to avoid the tt- and -ftT-decay bounds. In both the 
supersymmetric and nonsupersymmetric cases we are led to the prediction N e g ~ 2 for 
primordial nucleosynthesis. Finally the gravitino in the supersymmetric model has mass 
< 1 keV and hence satifies cosmo logical constraints. 

6. Naturalness 

We have seen that both the oscillation-decay solution for both the solar- and atmo- 
spheric neutrino problems, and the scalar-emission explanation of the excess double beta 
decay events, point to the necessity for scalar v. e.v.'s of about 1 keV. Generically a small 
scalar mass and expectation value are unacceptable because they are not stable under 
renormalization down from higher scales. We wish, in this section, to quantify this state- 
ment and to determine what constraints are required for the parameters of the theory in 
order to acheive this stability, or naturalness. 

The small parameters which are required to ensure such small scalar masses and v . e.v .'s 
appear in the scalar potential, which in the present model has the following renormalizable 
form: 

V{cj> e , eM , H) = -^M 2 ~ fil\^\ 2 - H 2 H rfH + U<Pe\ 4 + + ^e\ 2 \^\ 2 

+ £ e „|0 e | 2 #t# + £ Mff |0 eM | 2 #t# + Chh(H^H) 2 . 

(22) 

The hierarchy problem arises because, whereas [i H ~ 100 GeV, /j e and /U M can only be 
~ 1 keV. Similarly, although this hierarchy does not preclude some of the dimensionless 
couplings, namely £ HH , £ ee , £e/x an d £ma« f rom being 0(1), those which couple the light to 
the heavy fields, ^ eH and cannot be larger than < 0{pi 2 / [i 2 H ) ~ 10~ 16 . Otherwise 
quartic interactions such as \4> e \ 2 H^H would generate large mass terms for (p e once H 
receives its v. e.v. of 174 GeV. 

We must ask whether such small values for the couplings are ruined when they are run 
down from higher scales. For definiteness we define our running within the 'decoupling- 
subtraction' renormalization scheme, using dimensional regularization [22] . In this scheme 
all couplings are run between particle threshholds using modified minimal subtraction, and 
each particle is integrated out as the renormalization point is lowered below its threshhold. 
Within this framework, the MS running between threshholds only introduces a logarithmic 
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dependence on heavy mass scales, and the potentially dangerous powers of heavy scales 
arise when these heavy particles are integrated out. 

There are several kinds of graphs to consider, depending on which particles circulate 
in the loop. The most dangerous ones are those which involve the heaviest particles that 
are coupled the most strongly to the light scalars. 

• Higgs Loops: It is fairly easy to see that the running of the small couplings due to loops 
involving the standard-model Higgs are just as small as are the couplings themselves. For 
example the graphs of Figs. (1) give the following contributions to 5/j, 2 and 5£i H : 

(<K 2 ) Fig . i ~ %^ ~ (0-1 keV) 2 , and (<^) Fig . 1 ~ ~ 10" 18 , (23) 

where we take M H ~ 100 GeV, and we ignore all logarithms of large mass ratios. 

• Neutrino Loops: The really dangerous graphs involve heavy neutrinos, since the couplings 
of these particles to the light scalars is constrained by the requirement that a sufficiently 
large v e — mass matrix be generated at low energies. For example, the graphs of Fig. (2) 
produce the following contribution to 5^1 and 5£ eH : 

h 2 M 2 h 2 a 2 
(^)pig. 2 ~ -j^r, and (<^) Fig . 2 ~ T £|, (24) 

with similar contributions to S/j,^ and S^h- 

We see here the difficulty with having neutrinos in the MeV mass range. We work with 
the neutrino masses of most interest for double-beta decay and nucleosynthesis, M Se ~ 350 
MeV and M Sm ~ 1.5 MeV. We've taken M Se > 350 to avoid all bounds from K -> ev and 
tv — > ev decays [14]. The mixing of such a massive neutrino with v e is then constrained 
principally from weak universality, and satsifies |V e sJ 2 < 0.05 [15]. The bounds we obtain 
are listed in Table II. From this table we see that even if all couplings are made as large 
as they can be while remaining consistent with naturalness (and direct phenomenological 
constraints), their contributions to v e — mass matrix still fall two orders of magnitude 
short of what is required. 

The naturalness constraint can be eased by allowing the neutrinos to become lighter, 
since this allows their couplings to become weaker while keeping the desired low-energy 
v e — mass matrix fixed. This can only be done at the expense of reintroducing the 
cosmological problems and giving up our explanation for the double-beta decay anomaly. 
We consider a different approach in the next section. 
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Source 



Constraint 



Weak Universality g e < 4.5 x 10 4 771350 

ix —> [iv g^< 8.6 x 10~ 7 mi. 5 

Naturalness (5/j%) h e , < 4 x 10~ 5 to^q 

Naturalness 9eh e , 9eh^, g^h^ < 1 x 10" 7 
z/ e - 1^ Mass Matrix #g/i e ~4x 10~ 10 

#e£ M fy" - 2 X 10 ~ 12 m 350 mi.5 



Table II 

A summary of the constraints on the scalar-neutrino couplings. 771350 denotes (M Se /350 MeV) and mi. 5 
represents (M Sfl /l.5 MeV) . 

7. A Supersymmetric Model 

An alternative approach is to hold the line on the heavy-neutrino masses, and to 
instead try to work within a framework for which the naturalness conditions take a weaker 
form. There are two ways known to make a very light scalar v.e.v. natural. One either 
demands that the light scalar be a composite of underlying nonscalar particles, or one 
arranges for cancellations between fermions and bosons in the renormalization of small 
parameters. These cancellations can be ensured naturally if the model is supersymmetric. 

We have tried both approaches, but have only found a workable example in the su- 
persymmetric class. Here the contributions of the superpartners of each of the particles 
introduced so far largely cancel in their contributions to 5£i H and This cancellation 
is only accurate enough for scalar particles at the keV scale if the relevant suprmultiplets 
are themselves split in masses only by roughly a keV. Yet this must be made consistent 
with the fact that the superpartners of the charged leptons and quarks must have masses 
that are of order the electroweak scale. 

This may be consistently done within a supersymmetric framework if these highly 
degenerate multiplets are much more weakly coupled to the supersymmetry-breaking sector 
than are the other garden-variety particles we know. We demonstrate that this is possible 
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in this section by presenting an existence proof, in the form of a supersymmetric extension 
of our renormalizable model of Section (4). We shall find that, within this extension, the 
naturalness requirements are significantly relaxed. 



7.1) The Model 

The model we shall consider consists of the straightforward supersymmetric extension 
of our renormalizable model of Section (4). We imagine, therefore, supplementing the 
supersymmetric standard model (SSM) with four electroweak-singlet left-chiral superfields, 
carrying the following U(l) e xU(l) At charge assignments: 



$e 3 {<Pe, Xe} ~ (~2, 0), $ eM 3 {<p efl , Xe M } ~ (-1, -1), 

Nt 3 {M?,8f} ~ (±1,0), AT± 3 {A/* a±} ~ (0,±1). 



(25) 



Our notation here lists the scalar (cf) e ) and spinor (% e ) components of each superfield (<& e )- 
The most general renormalizable U(l) e xU(l) M -invariant superpotential becomes: 

W = W SSM + g e V e N-H u + g^»N-H u + ^-N+N+$ e + h fI N+N+^ efl 
+ m Se N+N- +m s ^N+N-, 

where 3 {Si, ipi} represent the SSM left-handed lepton-doublet superfields, which trans- 
form under U(l) e xU(l) M as * e ~ (1,0) and ~ (0, 1). H u similarly represents the SSM 
Higgs superfield whose v.e.v. gives up-type quarks their masses. 

This model as it stands does not break supersymmetry, and in the limit of negligible 
lepton- number breaking v.e.v.'s for the scalar fields, e , (p efI: Nf- and Af^ : it predicts com- 
pletely degenerate supermultiplets having the following masses: Four massless multiplets 
($ e , $ eAt , and one linear combination of and the neutrino multiplet in for both 
i = e and fx); a pair of degenerate massive multiplets (N~ and the other combination of 
iV+ and the neutrino multiplet in \l/ e ) having masses M Se = ^m 2 Se + g 2 v 2 ; and a similar 
pair for the muon-neutrino sector. Notice that, because of the unbroken supersymmetry, 
the scalars <j) e and e/i are massless, regardless of how large the couplings g e and h e should 
be. 

All of the interactions that are required to account for the solar- and atmospheric- 
neutrino problems, as well as the double-beta decay excess, appear in the superpotential of 
eq. (26) once this model is written out in terms of its component fields. This is because this 
superpotential contains a counterpart for each of the terms in the renormalizable couplings 
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of eq. (12), except for those terms which involve 0* or 4>* e ^ coupled to left-handed fields. 
Since these last terms do not appreciably affect the v e — sector, this difference is not 
important for our purposes. 



7.2) Supersymmetry Breaking 



In order to proceed further, some properties must be assumed for the supersymmetry- 
breaking sector of the model. Most importantly, it is necessary to determine whether it 
is possible to break supersymmetry in such a way as to give electroweak- scale masses to 
the superpartners of the ordinary quarks, leptons and gauge fields, and yet to still keep 
the mass-splittings within the supermultiplets $ e , <E> eM , and sufficiently small to 
permit naturally small expectation values (0 e ) ~ {4>e^) ~ 1 keV. 

In order to address this question, we imagine that supersymmetry is spontaneously 
broken when the auxiliary field of some electroweak- and U(l) e xU(l) M singlet superfield, 
Y, develops an expectation value: (Y) = A^#^, where 9 L is the Grassman coordinate of 
Y. Since (Y) 2 = 0, the most general nonderivative couplings this expectation value may 
have with the other fields are given by replacing Y — > (Y) in the following Kahler and 
superpotentials: 



SK = A(X,X*)Y*Y +[B(X,X*)Y* + c.c], and SW = C(X) Y. (27) 



In these expressions X a 3 {X a , x a } collectively represent all of the superfields of the SSM 
as well as the additional supermultiplets that were introduced in eq. (25) above. Once Y 
is replaced by (Y), 5K and 5W reduce to a set of soft SUSY- breaking contributions to 
the scalar potential, SV, as well as to the Yukawa-couplings, 5£ yu k- If F£ represents the 
auxiliary field for supermulitplet X a , then these contributions may be written 



5V(X,X*) =A 4 S A{X,X*) +A 



Q R 

_F» + C(*) + c.c. 



_A^ ( d 2 B 
yuk 2 \dX-dX b 



(28) 



X a lLX b + (gaugino terms) + c.c. 



We adopt the usual 'hidden-sector' scenario [23] in which the SUSY-breaking field, 
Y, couples to electroweak-scale fields only through interactions which are suppressed by a 
large mass scale, M. We may then write the most general lowest-dimension contributions 
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which involve the non-SSM fields as: 2 



6A 



5B 



M 2 
1 

M 



a e |$ e | 2 + a 



MI $ eM | 2 +«iV+iV- 



+ a'^N+N- +c.c 



)+££«? 



i=e,/j, 



n=± 



6 R |$ e | 2 



+ 6 M |$ e ^| 2 + (6'eiV+iV- + ft^+JV" + c.c.) + £ £ 61 



I AT, 



n\2 



5C = c e N+N~ + c^N+N- 



1 

+ M 



e/i 



(29) 

Using these expressions in eqs. (28) gives the induced SUSY-breaking interactions for our 
model. The induced scalar mass terms which arise in this way are generically 0(A 2 /M), 
with the exception of the terms c e N+N~ + c^N+N~ in 5C, which induce masses for 
that are 0(A S ). Since this latter size is much too large — as is the corresponding SSM 
term oc H u Hd — we take c e = = 0. Since SC is a contribution to the superpotential, 
this choice is technically natural due to the celebrated perturbative nonrenormalization 
theorems of supersymmetry. 

Generally the heavy mass scale, M, is chosen to ensure that the generic scalar masses 
that are obtained for the superpartners of the quarks and leptons are of order the weak 
scale, A 2 /M ~ v. With this choice the induced masses for scalars like 4> e would also be 
0(v), which is unacceptably large. We must therefore suppose that all contributions to 
8A through 8C are much smaller than would be indicated simply by the powers of M in 
eqs. (29) . We choose to parameterize this extra suppression by taking M to be much larger 
than the heavy mass, M SSM , which appears in the corresponding SSM terms (which we 
B SSM and C SSM ). This might be arranged in an underlying theory by having all 



call A 



SSM 5 



of the couplings between the supersymmetry-breaking sector and the very light fields, $j 
and N^, involve the virtual exchange of particles of mass M , while the coupling between 
Y and the SSM fields arises through the exchange of lighter particles of mass M SSM . M, 
M SSM and A s must therefore be related by: A 2 /M ~ 1 keV and A 2 /M SSM ~ 100 GeV, so 



M/M s 
then M 



SSM 



(100 GeV/1 keV) 
~ 10 11 GeV and A, 



10 8 . If we choose, for example, M 
t 3 x 10 6 GeV. 



10 19 GeV, 



The million-dollar question is whether this hierarchy of SUSY-breaking terms is tech- 
nically natural. The potentially difficult hierarchy to maintain is that between the sizes of 



If M should be as large as the Planck scale, then there would also be additional terms which arise 
from the elimination of the auxiliary fields of the supergravity multiplet. 



18 



8 A and A SSM , and of SB and B SSM since these are part of the Kahler potential, and so are 
not protected by the nonrenormalization theorems. This hierarchy would certainly be sta- 
ble if all of the couplings between SSM particles and the very light sector were suppressed 
by factors like M SSM /M, but this is not true for our model due to the comparatively large 
couplings, <7j and hi, in the supersymmetric superpotential of eq. (26). We must therefore 
check whether these interactions can communicate the large SSM SUSY-breaking scale to 
the light supermultiplets. 

Some of the potentially dangerous graphs are given in Fig. (3). We write these graphs 
using the supergraph formulation, in terms of which the supersymmetric cancellations are 
the most explicit. Fig. (3a) relates the splitting in the Nf 1 supermultiplets to those in $ e 
and $ eM , and gives: 



from which we infer: 5^, ~ (h 2 /16n 2 ) (a+A^/M 2 ). This has the same form as for the 
nonsupersymmetric model, except M 2 e has been replaced by the splitting within the iV+ 
supermultiplet: AM 2 ~ a+ K A S /M 2 . This is naturally small enough for our purposes, since 
the direct contribution fi 2 ~ a e A^/M 2 is by assumption 0(1 keV). But it underlines the 
fact that the members of the Nf 1 supermultiplets cannot be split in mass by more than a 
few orders of magnitude more than are the very light states in $ e or $ eAt . 

Fig. (3b) illustrates the other potentially dangerous combination, in which it is the 
mass splitting, Mj = A^(a SSM /M 2 SM ), between the SSM sneutrinos and neutrinos that 
sets the scale of the important loop momenta. This graph produces the following estimate 
for 5a e : 



from which we obtain: b~[i\ ~ (/i 2 /167r 2 ) (g e v/m s J 2 (m Se /Ms) 4 Mj This has a simple 
interpretation in terms of ordinary graphs involving the component fields. It reflects the 
fact that the contribution of the massive sneutrino is proportional to (/i 2 ff /167r 2 ) Mj, 
where h e s ~ h(g e v /m s J(m 2 e /Mj) is the effective coupling between the heavy sneutrino 
and the light scalar. The factor of g e v/m Se reflects the mixing angle between the SSM 
sneutrino multiplet, \I/ e , and the multiplets, A^, which couple to e . The remaining factor, 
m 2 e /Mj, similarly reflects the SUSY-breaking mixing angle between the heavy sneutrino 
mass eigenstate and the lighter scalar flavour eigenstates. 

We see, then, that this is indeed a stable hierarchy of scales, because the loop-induced 
splitting in the light multiplets due to the weak-scale splitting of the SSM multiplets is 




(30) 




(31) 
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sufficiently suppressed by the small coupling between these two sectors. 

8. Conclusions 

We have constructed several models which incorporate the phenomenological oscilla- 
tion/decay solution of Acker et.al. [5] for the solar and atmospheric neutrino problems. In 
so doing we have shown how the neutrino-scalar couplings that are required to produce the 
desired rate for neutrino decay, are also of the strength that is required to account for the 
recently-observed excess double beta decay events. This excess would then be interpreted 
as the emission of light scalars in addition to the electrons, with a predicted sum electron 
spectrum that falls into the 'Ordinary-Majoron' class of Ref. [9]. We regard this potential 
integration of the double-beta decay anomaly with the mechanism underlying the solar 
and atmospheric problems to be of sufficient interest to justify putting aside the poor 
agreement with the combined Homestake and GALLEX data, and to motivate a closer 
inspection of the implications of this class of models. 

In order to investigate this framework more completely, it must first be embedded into 
a renormalizable model that can be applied up to the weak scale. It is only then that it 
may be confronted with the complete range of experimental information that is available 
to constrain potential neutrino physics. We have made this comparison and find that the 
most stringent such constraints come from three sources: (a) Nucleosynthesis, (b) Double 
Beta Decay, and (c) Naturalness. We address each of these in turn. (They are summarized 
in Table III.) 

Any model for which neutrino decays play a role in explaining the solar-neutrino 
deficit almost inevitably involves very light degrees of freedom which couple appreciably 
to the electron neutrino. They are therefore typically hard pressed to account for the 
success of the standard Big-Bang framework of primordial nucleosynthesis, which tightly 
constrains the number of relativistic degrees of freedom that can be present when the 
temperature of the universe is of order (0.1 - 1) MeV. A conservative attitude is to simply 
ignore this constraint, with the rationale that the uncertainties involved in predicting the 
behaviour of the universe at such an early epoch are potentially quite large. Although we 
are basically sympathetic to this point of view, we have chosen here to see what is required 
to accomodate this remarkable success of standard cosmology. We find that consistency is 
possible, along the lines of the mechanism suggested in Ref. [19], provided that a relatively 
heavy state (such as a sterile neutrino) exists with a mass around 1.5 MeV. In this case 
we predict iV efr ~ 2, in good agreement with the data. 

The double-beta decay excess events also point to a heavy neutrino state, although in 
this case potentially much heavier than a few MeV. This condition arises as a special case 
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Criteron 



Constraint 



Atmospheric, solar-z/ anomalies 
Double beta decay 
v e , mixing with sterile v 
Cosmology 
Naturalness 



Table III 

A summary of the constraints which must be satisfied by the class of models defined in section 2 and the 
renormalizable model of section 4 in particular. The lepton number carrying scalars 4> e ,4>e^ are defined 
in section 2. The lepton number carrying heavy sterile neutrinos S e ,S M and their masses M Se ,M Sfl are 
defined in section 4. Only the supersymmetric generalization of the model in section 4 can incorporate all 
the criterea listed in column one and it does so by removing the naturalness constraint as given in the last 
row above. 

of a more general result which holds for any renormalizable model for which all light scalars 
are electroweak singlets [9]. The result holds because for these types of models no direct 
coupling is possible between the electron neutrino and the light scalar that is emitted in 
the decay. The act of emission therefore only occurs because of the mixing that arises due 
to the mismatch between the neutrino electroweak and mass eigenstates. The decay rate 
therefore goes to zero in the limit that all neutrinos are degenerate, and so in particular 
when all neutrinos are light in comparison with the MeV-scale energies that characterize 
the decay. 

By contrast with the previous considerations which point toward the existence of heavy 
neutrinos with masses in the MeV range, naturalness considerations generically prefer all 
new particle states to be lighter than this scale. Otherwise loops in which these heavy 
particles circulate produce contributions to the light scalar potential that are much larger 
than those that are permitted if the scalars are to have v.e.v.'s that are 0(1 keV). In the 
generic case, as may be seen from Table II and III, this naturalness criterion conflicts with 
the requirements of nucleosynthesis and double-beta decay, although only by one or two 



(0 e ) , Wen) ~ 1 keV 
(0 e ) , (0 eM ) ~ 1 keV , M Se > 60 MeV 
M Se < 10 MeV or M Se > 350 MeV 
M Se > 2.3 MeV , 0.7 MeV < M Sm < 2.3 MeV 
M Se ,M SM <lMeV 
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orders of magnitude. 

If these naturalness considerations are eschewed, our renormalizable model of Section 
(4) satisfactorily accounts for the solar and atmospheric neutrino deficits, explains the 
excess double-beta decay events, and satisfies all phenomenological constraints, including 
those from primordial nucleosynthesis. It robustly predicts that future solar-neutrino ex- 
periments will find some oscillations from electron- to muon-type neutrinos, together with 
an overall depletion of the solar flux due to neutrino decays with a lab-frame lifetime in 
the vicinity of 1000 sec. It also predicts that lower-energy neutrinos are more depleted 
than are higher-energy ones, requiring either the Homestake or the GALLEX signals to 
presently be incorrect at the several a level. 

In the particular renormalizable realization that is explored here, the model also pre- 
dicts the existence of sterile neutrinos in both the few MeV and the several-hundred MeV 
ranges. The MeV neutrino is appreciably mixed with z/ M , and the hundred MeV neutrino 
mixes significantly with v e . Both types of mixings may be amenable to detection through 
closer scrutiny of the 7r — > [iv decay spectrum, and in weak-universality measurements of 
the electron charged-current interaction strength. 

Taking the naturalness requirement seriously, one may simply accept very light neu- 
trino states, by ignoring the conflict with nucleosynthesis and the double-beta decay 
anomaly. Alternatively, if the low-energy model is embedded within a supersymmetric 
framework, all three types of constraints may be satisfied. In particular, we have shown 
how supersymmetry can be consistent with weak-scale masses for the superpartners of all 
presently-known particles, as well as keV-scale splittings among the new light supermulti- 
plets. This keV-scale splitting is what sets the size of the light scalar potential, and ensures 
the proper coupling strength for these scalars in the low-energy phenomenological model. 

In any of these scenarios we find that the atmospheric, solar and double beta decay 
anomalies lead us towards new scalars at very low energy scales, and new neutrino physics 
at MeV scales. 
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9. Figure Captions 



• Figure (1): One-loop contributions to the small scalar couplings which arise due to 
couplings with the heavy Standard-Model scalar. 

• Figure (2): Dangerous one- loop contributions to the small scalar couplings due to loops 
with heavy neutrinos. 

• Figure (3): The supergraphs which most strongly couple the SUSY-breaking scale of the 
SSM to the light multiplets. 
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